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Heart  rate  variability  (HRV)  decreases  during  hemorrhage,  and  has  been  proposed  as  a 
new  vital  sign  to  assess  cardiovascular  stability  in  trauma  patients.  The  purpose  of  this 
study  was  to  determine  if  any  of  the  HRV  metrics  could  accurately  distinguish  between 
individuals  with  different  tolerance  to  simulated  hemorrhage.  Specifically,  we  hypothesized 
that  (1)  HRV  would  be  similar  in  low  tolerant  (LT)  and  high  tolerant  (HT)  subjects  at  pre¬ 
syncope  when  both  groups  are  on  the  verge  of  hemodynamic  collapse;  and  (2)  HRV  could 
distinguish  LT  subjects  at  presyncope  from  hemodynamically  stable  HT  subjects  (i.e.,  at  a 
submaximal  level  of  hypovolemia).  Lower  body  negative  pressure  (LBNP)  was  used  as  a 
model  of  hemorrhage  in  healthy  human  subjects,  eliciting  central  hypovolemia  to  the  point 
of  presyncopal  symptoms  (onset  of  hemodynamic  collapse).  Subjects  were  classified  as 
LT  if  presyncopal  symptoms  occurred  during  the  —15  to  —60  mmHg  levels  of  LBNR  and  HT 
if  symptoms  occurred  after  LBNP  of  —60  mmHg.  A  total  of  20  HRV  metrics  were  derived 
from  R-R  interval  measurements  at  the  time  of  presyncope,  and  at  one  level  prior  to  pre¬ 
syncope  (submax)  in  LT  and  HT  groups.  Only  four  HRV  metrics  (Long-range  Detrended 
Fluctuation  Analysis,  Forbidden  Words,  Poincare  Plot  Descriptor  Ratio,  and  Fractal  Dimen¬ 
sions  by  Curve  Length)  supported  both  hypotheses. These  four  HRV  metrics  were  evaluated 
further  for  their  ability  to  identify  individual  LT  subjects  at  presyncope  when  compared  to 
HT  subjects  at  submax.  Variability  in  individual  LT  and  HT  responses  was  so  high  that  LT 
responses  overlapped  with  HT  responses  by  85-97%. The  sensitivity  of  these  HRV  metrics 
to  distinguish  between  individual  LT  from  HT  subjects  was  6-33%,  and  positive  predictive 
values  were  40-73%. These  results  indicate  that  while  a  small  number  of  HRV  metrics  can 
accurately  distinguish  between  LT  and  HT  subjects  using  group  mean  data,  individual  HRV 
values  are  poor  indicators  of  tolerance  to  hypovolemia. 
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INTRODUCTION 

In  trauma  patients,  the  primary  cause  of  death  within  the  first 
hour  of  injury  is  hemorrhage  (Bellamy,  1984;  Sauaia  et  al.,  1995; 
Champion  et  al.,  2003).  Early  detection  of  hypovolemia  in  trauma 
patients  is  therefore  critical  to  improving  triage  and  providing 
life  saving  interventions  (LSIs).  Unfortunately,  the  severity  of 
hemorrhage  or  central  hypovolemia  is  difficult  to  detect  by  first 
responders  or  other  emergency  medical  personnel,  as  the  tra¬ 
ditional  vital  signs  currently  available  to  them,  such  as  arterial 
pressure,  heart  rate,  respiration  rate,  and  pulse  character  do  not 
consistently  change  until  hemorrhage  has  progressed  to  the  point 
of  decompensation  (i.e.,  hemodynamic  collapse;  Heckbert  et  al., 
1998;  Parks  et  al.,  2006;  Convertino  et  al.,  2008).  Assessing  the 
severity  of  hemorrhage  in  trauma  patients  is  further  complicated 
by  the  reality  that  there  are  individuals  who  have  a  low  tolerance 
to  hemorrhage  such  that,  for  the  same  level  of  blood  loss,  low  tol¬ 
erant  (LT)  individuals  reach  the  point  of  cardiovascular  collapse 
in  less  time  than  those  who  are  high  tolerant  (HT;  Shoemaker 


et  al.,  1973,  2006;  Klemcke  et  al.,  2011).  Thus,  there  is  an  urgent 
need  to  identify  physiological  parameters  that  can  accurately  track 
hypovolemia  in  trauma  patients,  and  have  the  sensitivity  to  pro¬ 
vide  early  identification  of  a  patient  that  has  a  low  tolerance  to  this 
physiological  insult.  New  vital  signs  are  needed  which  can  detect 
the  onset,  progression,  and  severity  of  hypovolemia  in  individual 
patients. 

The  measurement  of  the  variation  in  R-R  intervals,  i.e.,  heart 
rate  variability  (HRV;  Malik,  1996),  can  provide  information 
regarding  overall  cardiovascular  status.  For  instance,  healthy  indi¬ 
viduals  with  normal  autonomic  function  have  high  HRV,  while 
individuals  with  a  stressed  cardiovascular  system  (via  activity,  dis¬ 
ease,  or  dysfunction)  can  be  identified  by  reduced  HRV  (Thayer 
and  Sternberg,  2006) .  HRV  has  been  extensively  studied  in  hospital 
settings  such  as  intensive  care  units  to  assess  cardiovascular  sta¬ 
tus  and  predict  mortality  in  trauma  patients  (Winchell  and  Hoyt, 
1996,  1997;  Grogan  et  al.,  2005;  Norris  et  al.,  2005,  2008;  Mor¬ 
ris  et  al.,  2006).  Considering  that  monitoring  electrocardiogram 
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(ECG)  signals  is  non-invasive,  can  be  accomplished  in  all  eche¬ 
lons  of  care  (including  the  pre-hospital  environment),  and  can 
potentially  provide  information  regarding  cardiovascular  status, 
the  utility  of  HRV  as  a  pre-hospital  triage  assist  tool  is  the  focus 
of  many  investigations  (Cooke  et  al.,  2006a, b;  Batchinsky  et  al., 
2007a;  Cancio  et  al.,  2008;  Ryan  et  al.,  2011). 

One  criterion  for  the  development  of  accurate  and  timely 
assessment  of  the  clinical  status  of  patients  with  hemorrhage  is 
the  capability  to  investigate  the  continuous  physiological  response 
to  progressive  reductions  in  central  blood  volume.  This  capability 
has  become  available  with  the  emergence  of  lower  body  negative 
pressure  (LBNP)  as  a  human  experimental  model  for  hemorrhage 
(Cooke  et  al.,  2004).  Cardiovascular  responses  during  LBNP  are 
reproducible,  and  are  consistent  with  those  observed  during  hem¬ 
orrhage  (Convertino  and  Sather,  2000;  Convertino,  2001;  Cooke 
et  al.,  2004).  In  addition,  LBNP  allows  for  the  assessment  of 
cardiovascular  responses  which  are  specific  to  hypovolemia  as 
confounding  factors  associated  with  trauma  such  as  pain  and 
inflammation  are  absent.  Thus,  the  controlled  experimental  envi¬ 
ronment  of  LBNP  allows  for  optimal  conditions  in  which  to  study 
HRV.  Using  the  LBNP  model  to  induce  a  clinical  endpoint  (i.e., 
presyncope),  human  subjects  can  be  identified  as  having  HT  or 
LT  to  central  hypovolemia  (Sather  et  al.,  1986;  Convertino  and 
Sather,  2000;  Rickards  et  al.,  2011).  If  HRV  metrics  reflect  the 
status  of  cardiovascular  stability  during  progressive  central  hypo¬ 
volemia,  then  LT  subjects  would  be  expected  to  display  greater 
reductions  in  HRV  compared  to  HT  subjects.  We  therefore  used 
LBNP  tolerance  as  a  model  to  test  the  hypotheses  that:  ( 1 )  HRV  at 
presyncope  is  similar  in  LT  and  HT  subjects  when  both  groups  are 
on  the  verge  of  cardiovascular  collapse,  and  (2)  LT  subjects  have 
lower  HRV  at  presyncope  compared  to  HT  subjects  who  remain 
hemodynamically  stable. 

MATERIALS  AND  METHODS 

SUBJECTS 

This  study  was  conducted  under  a  protocol  reviewed  and  approved 
by  the  Brooke  Army  Medical  Center  and  US  Army  Medical 
Research  and  Materiel  Command  Institutional  Review  Boards  and 
in  accordance  with  the  approved  protocol.  All  studies  were  con¬ 
ducted  at  the  U.S.  Army  Institute  of  Surgical  Research,  Fort  Sam 
Houston,  TX,  USA.  Normotensive,  non-smoking  healthy  human 
volunteers  (N  =  120;  male  =  75,  female  =  45,  age  28  ±  1  years; 
height  174  ±  lcm;  weight  76  ±  lkg)  participated  in  this  study 
after  evaluation  of  their  medical  history  and  physical  examina¬ 
tion  by  a  physician  to  ensure  the  absence  of  previous  and  current 
medical  conditions  that  would  exempt  them  as  participants.  All 
female  subjects  were  confirmed  as  not  pregnant  by  a  urine  preg¬ 
nancy  test  conducted  within  an  hour  before  the  experiment.  For 
24  h  prior  to  the  study,  participants  were  instructed  to  maintain 
their  normal  sleep  patterns,  refrain  from  exercise,  and  abstain 
from  alcohol  and  autonomic  stimulants  including  prescription 
and  non-prescription  drugs  such  as  caffeine,  alcohol,  and  decon¬ 
gestants.  The  potential  effect  of  caffeine  withdrawal  on  baseline 
hemodynamics  was  not  assessed  in  this  study.  Subjects  were  given 
a  written  description  of  the  experimental  protocol  and  the  risks 
associated  with  the  study.  On  the  day  of  the  study,  the  subjects 
were  made  familiar  with  the  laboratory,  given  a  verbal  briefing  on 


the  protocol  and  procedures,  and  encouraged  to  ask  questions  of 
the  investigators.  All  subjects  signed  an  IRB  approved  informed 
consent  form. 

LBNP PROTOCOL 

Subjects  were  instrumented  for  standard  lead  II  ECG  to  record 
R-R  intervals  (RRI),  and  a  finger  cuff  to  record  beat-by-beat  fin¬ 
ger  arterial  pressure  by  photoplethysmography  (Finometer  Blood 
Pressure  Monitor,  TNO-TPD  Biomedical  Instrumentation,  Ams¬ 
terdam,  The  Netherlands).  To  simulate  hemorrhage  in  conscious 
humans,  central  hypovolemia  was  induced  by  application  of  LBNP. 
Previous  studies  have  shown  that  the  LBNP  protocol  can  closely 
simulate  the  hemodynamic  challenges  associated  with  pre-shock 
hemorrhage  (Cooke  et  al.,  2004;  Ward  et  al.,  2010).  Subjects  were 
positioned  supine  within  an  airtight  chamber  that  was  sealed  at 
the  level  of  the  iliac  crest  by  a  neoprene  skirt.  The  LBNP  protocol 
consisted  of  a  5-min  control  period  (baseline)  followed  by  5  min 
of  chamber  decompression  at  —15,  —30,  —45,  and  — 60mmHg, 
and  then  additional  increments  of  — 10  mmHg  every  5  min  until 
the  onset  of  cardiovascular  collapse  followed  by  a  10-min  recovery 
period.  The  application  of  LBNP  was  terminated  based  on  three 
criteria.:  (a)  sudden  onset  of  relative  bradycardia,  (b)  progres¬ 
sive  fall  of  systolic  pressure  below  80  mmHg;  and,  (c)  voluntary 
subject  termination  due  to  the  onset  of  presyncopal  symptoms 
such  as  sweating,  nausea,  dizziness,  vision  alterations,  or  general 
discomfort.  During  the  LBNP  protocol,  arterial  pressure  wave¬ 
form  (Finometer)  and  ECG  signals  were  sampled  continuously 
at  500  Hz  and  digitally  recorded  with  WinDaq  data  acquisition 
software  (DATAQ  Instruments,  Akron,  OH,  USA). 

HIGH  TOLERANCE  VS.  LOW  TOLERANCE 

Subjects  were  categorized  as  LT  or  HT  based  on  the  level  of  LBNP 
at  which  they  experienced  the  onset  of  cardiovascular  collapse 
(Sather  et  al.,  1986;  Rickards  et  al.,  2011).  LT  subjects  experienced 
cardiovascular  collapse  before  or  during  —60  mmHg  of  LBNP,  and 
HT  subjects  experienced  cardiovascular  collapse  after  —60  mmHg 
of  LBNP. 

DATA  ANALYSIS 

Arterial  pressure  waveform  and  ECG  signals  were  analyzed  during 
the  last  3  min  of  each  LBNP  level  using  a  commercially  available 
software  program  (WinCPRS,  Absolute  Aliens,  Turku,  Finland). 
Beat-to-beat  stroke  volume  was  derived  from  the  arterial  pressure 
waveform  using  the  pulse  contour  method  (Jansen  et  al.,  1990). 

Electrocardiogram  signals  were  visually  inspected  and  con¬ 
firmed  to  be  free  of  noise,  persistent  ectopic  beats,  or  arrhythmias. 
When  ectopic  beats  were  identified,  the  presence  of  a  single  event 
within  the  3 -min  analysis  window  resulted  in  the  correction  of 
the  beat  by  linear  interpolation  prior  to  analysis.  If  more  than 
one  ectopic  beat  occurred  during  an  LBNP  level,  the  ECG  data 
for  that  level  were  not  analyzed.  Acceptable  ECG  signals  were 
used  to  calculate  RRI  and  heart  rate.  The  RRI  signal  was  used 
to  derive  20  indices  of  HRV  described  in  detail  in  previous  pub¬ 
lications  from  this  laboratory  (Rickards  et  al.,  2010a;  Ryan  et  al., 
2010).  Linear  methods  were  used  to  calculate  HRV  metrics  in 
both  the  time  domain  [«  =  9;  RRI  standard  deviation  (RRISD), 
RRI  root  mean  squared  standard  deviation  (RMSSD),  percent¬ 
age  adjacent  RRIs  varying  by  at  least  50  ms  (pNN50),  Poincare 
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plot  descriptors  standard  deviation  1  and  standard  deviation  2 
(SD1  and  SD2),  SD1/SD2  ratio,  SD2/SD1  ratio,  complex  demod¬ 
ulation  high  frequency  (CDM  HF),  complex  demodulation-low 
frequency  (CDM  LF)  ] ;  and  frequency  domain  [  n  =  2;  RRI  low  fre¬ 
quency  power  (RRI  LF),  RRI  high  frequency  power  (RRI  FiF) ] . 
For  frequency  domain  analysis,  the  RRI  signals  were  plotted  using 
linear  interpolation,  sampled  at  5  Hz,  and  then  passed  through  a 
low-pass  impulse  response  filter  with  a  cutoff  frequency  of  0.5  Hz. 
Data  sets  were  submitted  to  a  Fourier  transform  with  a  Hanning 
window.  Non-linear  methods  were  used  to  calculate  heart  rate 
complexity  metrics  from  the  RRI  signal  (n  =  9;  sample  entropy 
(SampEn),  Lempel-Ziv  entropy  (LZEn),  fractal  dimensions  by 
curve  length  (FD-L),  fractal  dimensions  by  dispersion  analysis 
(FD-DA),  symbol  dynamics  entropy  (SymDyn),  normalized  sym¬ 
bol  dynamics  entropy  (DisnEn),  long-range  detrended  fluctuation 
analysis  (DFA  long),  short-range  detrended  fluctuation  analysis 
(DFA  short),  and  Forbidden  Words  (FW). 

Heart  rate  variability  metrics  were  compared  in  LT  and  HT 
subjects  at  the  level  of  LBNP  at  which  cardiovascular  collapse  was 
imminent  (presyncope).  In  HT  subjects,  HRV  metrics  were  also 
calculated  one  level  of  LBNP  prior  to  the  presyncopal  level  (sub¬ 
max),  and  compared  with  LT  presyncopal  values.  HRV  metrics 
which  measured  equal  variability  in  LT  and  HT  at  presyncope  and 
also  less  variability  in  LT  at  presyncope  compared  to  HT  at  sub¬ 
max  were  further  evaluated  for  their  ability  to  distinguish  between 
individual  LT  and  HT  subjects. 

SENSITIVITY  AND  POSITIVE  PREDICTIVE  VALUE 

The  effectiveness  of  HRV  metrics  to  accurately  determine  group 
membership  (HT  or  LT)  was  evaluated  by  calculating  sensitivity 
and  positive  predictive  value.  Sensitivity  assesses  the  ability  of  a 
HRV  metric  to  correctly  identify  a  LT  subject  from  all  the  sub¬ 
jects  that  were  actually  LT  (i.e.,  true  positives/true  positives  +  false 
negatives).  Positive  predictive  value  is  a  measure  of  the  ability  of 
a  HRV  metric  to  identify  a  subject  as  LT  from  all  the  subjects  that 
“appear”  to  be  LT  (i.e.,  true  positives/true  positives  +  false  posi¬ 
tives).  A  logistic  regression  model  was  used  to  evaluate  the  validity 
of  HRV  metrics  as  predictors  of  low  tolerance  to  hemorrhage. 
Using  this  model,  cutoff  values  were  determined  for  each  parame¬ 
ter  at  each  level  of  LBNP  based  on  receiver-operator  characteristic 
analysis  with  specificity  no  less  than  0.95.  Sensitivity  and  posi¬ 
tive  predictive  values  were  then  calculated  based  on  this  criterion 
(Pregibon,  1981). 

STATISTICAL  ANALYSIS 

Data  are  expressed  as  mean  ±  SE.  Statistical  analysis  was  con¬ 
ducted  with  commercially  available  software  (SigmaStat;  Systat 
Software,  Richmond,  CA,  USA).  Unpaired  f-test,  Chi-square,  and 
z-test  analyses  were  used  to  compare  subject  demographics  and 
conditions  at  LBNP  termination  between  tolerance  groups.  For 
hemodynamic  responses,  a  two  way  analysis  of  variance  with 
repeated  measures  on  one  variable  was  used  to  determine  signifi¬ 
cant  variable  (tolerance  and  LBNP  level)  main  effects  ( p  <  0.05). 
Subsequent  multiple  comparison  tests  (Holm-Sidak)  determined 
significant  differences  between  HT  and  LT  groups  (p  <  0.05).  HRV 
measurements  at  presyncope  and  submax  in  LT  and  HT  subjects 
were  compared  by  unpaired  f-test. 


RESULTS 

SUBJECTS 

The  LBNP  procedure  was  conducted  on  120  subjects.  LBNP  was 
terminated  in  2  subjects  during  —  30mmHg,  in  5  subjects  during 
—45  mmHg,  in  26  subjects  during  —60  mmHg,  in  43  subjects  dur¬ 
ing  —70  mmHg,  in  29  subjects  during  —80  mmHg,  in  13  subjects 
during  —90  mmHg,  and  in  2  subjects  during  —100  mmHg.  As  a 
result,  33  subjects  were  identified  as  LT  (27.5%)  and  87  subjects 
were  identified  as  HT  (72.5%).  The  demographics  and  baseline 
mean  arterial  pressure  and  heart  rates  of  these  two  groups  of  sub¬ 
jects  were  not  statistically  different  (Table  1).  However,  baseline 
RRI  was  significantly  lower  in  LT  subjects  compared  to  HT  sub¬ 
jects.  The  presyncopal  conditions  observed  at  LBNP  termination 
are  shown  in  Table  2.  Cardiovascular  conditions  include  systolic 
pressure  below  80  mmHg  and/or  sudden  onset  of  relative  brady¬ 
cardia,  and  presyncopal  symptoms  represent  discomfort  reported 
by  the  subject.  LT  and  HT  subjects  experienced  similar  presyncopal 
conditions  at  LBNP  termination. 

RRI  TRACKS  HYPOVOLEMIA 

Stroke  volume  and  RRI  at  baseline  (LBNP  =  0  mmHg)  and  dur¬ 
ing  LBNP  in  LT  and  HT  subjects  are  shown  in  Figures  1A,B. 
LT  subjects  experienced  presyncopal  conditions  and/or  symptoms 
before  or  during  LBNP  level  of  —60  mmHg.  HT  subjects  experi¬ 
enced  presyncopal  conditions  and/or  symptoms  after  —60  mmHg. 


Table  1  |  Demographics  and  baseline  hemodynamic  variables  for  low 
tolerant  and  high  tolerant  subjects. 


Low  tolerant  High  tolerant  P  value 

N 

33 

87 

<0.001 

Male/female 

18/15 

57/30 

0.37 

Age  (years) 

27  ±1 

29  ±  1 

0.32 

Height  (cm) 

173  ±2 

174±1 

0.60 

Weight  (kg) 

76  ±3 

76  ±  1 

0.98 

Mean  arterial  pressure  (mmHg) 

99  ±2 

97  ±  1 

0.23 

Heart  rate  (beats/min) 

69  ±2 

64  ±  1 

0.06 

R-R  interval  (ms) 

896  ±22 

973  ±17 

0.01 

Data  are  expressed  as  mean  ±  SE. 

Table  2  |  Prevalence  of  presyncopal  conditions  at  LBNP  termination  in 
low  tolerant  and  high  tolerant  subjects. 

Conditions  at  LBNP 

termination 

Low  tolerant 

(n  =  33) 

High  tolerant 
(n  =  87) 

P  value 

Cardiovascular±  symptoms 

17  (52%) 

51  (59%) 

0.622 

Cardiovascular  only 

5  (15%) 

12  (14%) 

0.917 

Symptoms  only 

11  (33%) 

24  (27%) 

0.694 

Cardiovascular  conditions  include  systolic  blood  pressure  below  80  mmHg  and/or 
sudden  bradycardia.  Symptoms  represent  presyncopal  symptoms  such  as  sweat¬ 
ing,  nausea,  dizziness,  vision  alterations,  or  general  discomfort  as  reported  by  the 
subject.  Data  are  expressed  as  number  of  occurrences;  group  percentages  are 
shown  in  parentheses. 
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LBNP  Level  (mm  Hg) 


FIGURE  1  |  Stroke  volume  (A)  and  R-R  interval  (B)  during  progressive 
lower  body  negative  pressure  (LBNP)  in  low  tolerant  subjects  (closed 
circles,  solid  line)  and  high  tolerant  subjects  (open  circles,  dashed  line). 

At  baseline,  there  are  33  low  tolerant  subjects  and  87  high  tolerant  subjects. 


LBNP  Level  (mm  Hg) 


The  number  of  subjects  at  the  subsequent  time  points  during  LBNP  decline 
as  subjects  experience  presyncopal  symptoms  at  various  levels  of  LBNP  Data 
are  shown  as  mean  ±  SE.  *  Denotes  group  differences  at  the  same  LBNP 
level  (p  <  0.05). 


Data  for  HT  subjects  are  shown  to  LBNP  level  of  —  80mmHg. 
Both  stroke  volume  and  RRI  decreased  in  a  similar  incremen¬ 
tal  fashion  during  the  progression  of  LBNP,  and  average  RRI 
measurements  reflect  and  are  correlated  with  the  hypovolemia 
induced  by  LBNP  (amalgamated  r2  =  0.988  between  stroke  vol¬ 
ume  and  RRI).  Stroke  volume  was  not  different  between  LT  and 
HT  groups,  but  RRI  was  significantly  lower  in  LT  subjects  than 
HT  subjects  at  rest  and  during  —15,  —30,  and  —45  mmHg  LBNP 
(P  <  0.01). 

HRV  MEASUREMENTS  AT  PRESYNCOPE  AND  SUBMAX  LEVELS  OF 
LBNP 

Figure  2  demonstrates  the  variation  in  RRI  in  a  single  LT  sub¬ 
ject  and  a  single  HT.  These  sample  tracings  of  RRI  were  obtained 
during  LBNP  level  of  —60  mmHg,  and  are  3  min  in  duration.  The 
variation  in  RRI  is  low  in  the  LT  subject  who  is  on  the  verge  of 
cardiovascular  collapse  (at  presyncope).  In  contrast,  the  variability 
in  RRI  is  high  in  the  HT  subject  who  is  hemodynamically  stable 
(at  submax). 

Heart  rate  variability  was  quantified  at  presyncope  in  all  LT 
and  HT  subjects,  and  at  submax  in  HT  subjects  by  20  HRV  met¬ 
rics.  These  results  are  shown  in  Table  3.  The  first  two  columns 
in  Table  3  indicate  the  parameter  or  HRV  metric  and  its  pro¬ 
portionality  to  variability.  All  HRV  metrics  are  quantitatively 
proportional  (P)  to  HRV  except  SD2/SD1,  DFA  long,  DFA  short, 
and  FW,  which  are  inversely  proportional  (I)  to  HRV.  The  val¬ 
ues  of  proportional  metrics  decrease  as  variability  decreases,  while 
the  values  of  inversely  proportional  metrics  increase  as  variabil¬ 
ity  decreases.  Columns  3,  4,  and  5  in  Table  3  display  the  data 
for  three  groups:  LT  at  presyncope,  HT  at  presyncope,  and  HT 
at  submax.  The  comparisons  of  variability  between  LT  at  presyn¬ 
cope  verses  HT  at  presyncope,  and  LT  at  presyncope  verses  HT  at 
submax  are  also  shown  in  Table  3,  in  columns  6  and  7,  respec¬ 
tively.  The  hypotheses  of  our  study  were  ( 1 )  that  at  presyncope 
HRV  would  be  equal  in  LT  and  HT  subjects  (column  6)  because 
both  groups  are  on  the  verge  of  cardiovascular  collapse;  and  (2) 


that  LT  subjects  would  have  less  HRV  at  presyncope  compared 
to  HT  subjects  who  remain  hemodynamically  stable  at  submax 
(column  7).  These  hypotheses  are  stated  in  Table  3  above  the 
comparisons  of  HRV  metrics  (i.e.,  “LT  =  HT,”  “LT  <  HT”).  When 
comparing  HRV  between  groups  it  is  important  to  keep  in  mind 
that  the  numerical  values  of  proportional  metrics  are  consistent 
with  variability  (i.e.,  high  value  indicates  high  variability),  but 
numerical  values  of  metrics  which  are  inversely  proportional  to 
variability  are  opposite  to  variability  (i.e.,  high  value  indicates  low 
variability).  At  presyncope,  HRV  was  equivalent  between  LT  and 
HT  groups  as  measured  by  7  HRV  metrics  (bold  type  in  column 
6  of  Table  3).  When  HRV  was  compared  between  LT  subjects 
at  presyncope  and  HT  subjects  at  submax,  8  metrics  were  less 
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Table  3  |  HR,  RRI,  and  HRV  metrics  in  LT  at  presyncope,  HT  at  presyncope,  and  HT  at  submax. 


Metric 

Proportional  (P), 

LT  at 

HT  at 

HT  at 

LT  @  Presyncope  vs. 

LT  @  Presyncope  vs. 

Inversely  proportional  (1) 

presyncope 

presyncope 

submax 

HT@  Presyncope 

HT  @  Submax 

ECG  SIGNAL 

HR 

97.3  ±3.0 

1 1 7.9  ±  2.3 

106.3  ±  2.1 

LT  <  HT 

LT  <  HT 

RRI 

643. 6±  22.0 

530. 6±  11.6 

587.1  ±113.4 

LT  >  HT 

LT  >  HT 

HRV  metrics 

Hypothesis 

Hypothesis 

LT  =  HT 

LT  <  HT 

TIME  DOMAIN 

RRISD 

P 

64.1  ±6.0 

38.9±  2.3 

36.7  ±  1.9 

LT  >  HT 

LT  >  HT 

RMSSD 

P 

23.9±  5.8 

9.6±0.9 

14. 8±  1.2 

LT  >  HT 

LT  >  HT 

pNN50 

P 

3.2  ±  1.1 

0.9±0.2 

2.5±0.6 

LT  >  HT 

LT  =  HT 

SD1 

P 

1 3.2  ±  1.9 

6.8±0.6 

10.5±0.9 

LT  >  HT 

LT  =  HT 

SD2 

P 

88.2  ±8.1 

54.0±3.2 

50.1  ±2.6 

LT  >  HT 

LT  >  HT 

SD1/SD2* 

P 

0.146±0.010 

0.129±0.010 

0.201  ±0.011 

LT  =  HT 

LT  <  HT 

SD2/SD1 

1 

7.80  ±0.50 

10.17±0.53 

5.85±0.24 

LT  >  HT 

LT  <  HT 

CDM  HF 

P 

12. 8±  1.7 

6.9±0.6 

11.7  ±  1.1 

LT  >  HT 

LT  =  HT 

CDM  LF 

P 

27.1  ±2.3 

20.1  ±1.7 

27.9  ±1.9 

LT  >  HT 

LT  =  HT 

FREQUENCY  DOMAIN 

RRI  HF 

P 

216.1  ±88.7 

55. 6±  10.6 

188.9  ±38.1 

LT  >  HT 

LT  =  HT 

RRI  LF 

P 

432.7  ±  91.4 

302.4±62.3 

742.7±  111.4 

LT  =  HT 

LT  =  HT 

COMPLEXITY  DOMAIN 

SampEn 

P 

0.707  ±0.048 

0.57±  0.029 

0.821  ±0.031 

LT  >  HT 

LT  =  HT 

LZEn 

P 

0.504±  0.026 

0.420±0.017 

0.566±0.016 

LT  >  HT 

LT  <  HT 

FD-L* 

P 

1.572  ±0.026 

1.529±  0.016 

1.701  ±0.011 

LT  =  HT 

LT  <  HT 

FD-DA 

P 

1.152±  0.020 

1.1 28  ±  0.01 2 

1. 1 86  ±  0.01 1 

LT  =  HT 

LT  =  HT 

SymDyn 

P 

0.521  ±0.018 

0.474  ±0.011 

0.580±0.010 

LT  >  HT 

LT  <  HT 

DisnEn 

P 

3.142  ±0.111 

2.846  ±0.067 

3.482  ±0.059 

LT  >  HT 

LT  <  HT 

DFA  long* 

1 

1.054  ±0.031 

1.100±0.025 

0.887  ±  0.022 

LT  =  HT 

LT  <  HT 

DFA  short 

1 

1.564  ±0.052 

1.619±0.034 

1.605±0.029 

LT  =  HT 

LT  =  HT 

FW* 

1 

67.3  ±0.9 

69.1  ±0.7 

63.8±0.7 

LT  =  HT 

LT  <  HT 

Data  are  shown  in  three  groups:  low  tolerant  (LT)  at  presyncope:  high  tolerant  (HT)  at  presyncope:  and  HT  at  submax.  Values  (mean±  SE)  are  shown  for  Heart  Rate 
(HR),  R-R  interval  (RRI),  RRI  standard  deviation  (RRISD),  RRI  root  mean  squared  standard  deviation  (RMSSD),  percentage  adjacent  RRIs  varying  by  at  least  50  ms 
(pNN50),  Poincare  plot  descriptors  standard  deviation  1  (SD1),  and  standard  deviation  2  (SD2),  SD1/SD2  ratio,  SD2/SD1  ratio,  complex  demodulation  high  frequency 
(CDM  HF),  complex  demodulation-low  frequency  (CDM  LF),  RRI  high  frequency  power  (HF),  RRI  low  frequency  power  (LF),  sample  entropy  (SampEn),  Lempel-Ziv 
entropy  (LZEn),  fractal  dimensions  by  curve  length  (FD-L),  fractal  dimensions  by  dispersion  analysis  (FD-DA),  symbol  dynamics  entropy  (SymDyn),  normalized  symbol 
dynamics  entropy  (DisnEn),  long-range  detrended  fluctuation  analysis  (DFA  long),  short-range  detrended  fluctuation  analysis  (DFA  short).  P  indicates  HRV  metrics 
values  proportional  to  RRI  variability:  I  indicates  HRV  metrics  values  inversely  proportional  to  RRI  variability.  Comparisons  between  groups  are  shown  as  LT=  HT 
LT  <  HT  and  LT>  HT:  <,  >  denote  significant  difference  between  LT  and  HT  group  means  (p<  0.05):  bold  type  comparisons  support  the  related  hypothesis.  *  Denotes 
the  metrics  with  responses  to  support  both  hypotheses. 


variable  in  LT  at  presyncope  (bold  type  in  column  7  of  Table  3). 
The  responses  of  DFA  long,  FW,  SD1/SD2,  and  FD-L  were  the 
only  metrics  to  support  both  hypotheses,  and  are  indicated  by 
an  asterisk  (*)  in  Table  3.  These  four  HRV  metrics  were  further 
evaluated  for  their  ability  to  distinguish  individual  LT  and  HT 
subjects. 

DFA  LONG.  FW.  SD1/SD2,  AND  FD-L  IN  INDIVIDUAL  LT  AND  HT 
SUBJECTS 

DFA  long,  FW,  SD1/SD2,  and  FD-L  in  individual  LT  subjects  at 
presyncope  and  individual  HT  subjects  at  submax  are  shown  as 
box  and  whisker  plots  in  Figures  3A-D.  Statistical  analysis  of 
group  averages  (Table  3)  reveal  significant  differences  in  the  four 


HRV  measurements  between  LT  and  HT  groups.  However,  when 
the  measurements  of  DFA  long,  FW,  SD1/SD2,  and  FD-L  from 
individual  subjects  are  graphically  displayed  (Figures  3A-D),  it 
is  evident  that  individual  responses  from  LT  subjects  at  presyn¬ 
cope  are  predominantly  within  the  range  of  responses  observed  in 
HT  subjects  at  submax  despite  the  fact  that  LT  subjects  are  on  the 
verge  of  cardiovascular  collapse  and  HT  subjects  are  hemodynam- 
ically  stable.  The  percentage  of  LT  responses  that  were  coincident 
with  HT  responses  were  97%  for  DFA  long  and  FW,  94%  for 
SD1/SD2,  and  85%  for  FD-L.  The  sensitivities  of  DFA  long,  FW, 
SD1/SD2,  and  FD-L  to  identify  individual  LT  subjects  were  low 
(12,  6,  24,  33%),  and  positive  predictive  values  were  50,  40,  67, 
and  73%. 
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FIGURE  3  |  (A)  Long-range  detrended  fluctuation  analysis  (DFA 
long),  (B)  forbidden  words  (FW),  (C)  Poincare  plot  standard 
deviations  ratio  (SD1/SD2),  and  (D)  fractal  dimensions  by  curve 
length  (FD-L)  in  low  tolerant  (LT)  at  presyncope  (solid  circles, 


n  =  33)  and  high  tolerant  (HT)  at  submax  (open  circles,  n=  87)  level 
of  lower  body  negative  pressure  (LBNP).  Data  are  shown  as  box 
(25th/75th  percentiles)  and  whisker  (90th/10th  percentiles)  plots 
with  median  value  (black  line). 


DISCUSSION 

To  investigate  the  application  of  HRV  monitoring  for  assessment 
of  hemodynamic  stability  during  hemorrhage,  we  used  an  exper¬ 
imental  model  to  compare  HRV  metrics  in  LT  and  HT  subjects 
during  LBNR  In  order  for  any  specific  HRV  metric  to  be  con¬ 
sidered  as  a  valid  candidate  for  the  assessment  of  cardiovascular 
stability  in  the  context  of  central  blood  volume  loss,  two  condi¬ 
tions  must  hold  true  using  this  model:  (1)  LT  and  HT  subjects 
should  display  similar  HRV  at  presyncope;  and  (2)  HRV  in  LT 
subjects  should  be  less  (i.e.,  LT  subjects  would  be  less  stable)  than 
HRV  in  HT  subjects  at  the  point  in  time  when  LT  subjects  expe¬ 
rienced  presyncope.  These  conditions  are  based  on  the  premise 
that  hemodynamic  instability  reflects  a  specific  physiological  con¬ 
dition  defined  by  the  inability  of  cardiovascular  mechanisms  to 
adequately  compensate  for  reduced  central  blood  volume  in  all 
subjects  independent  of  their  tolerance,  and  that  LT  subjects 
are  hemodynamically  unstable  at  presyncope  when  HT  subjects 
remain  stable.  Our  results  indicate  that  of  the  20  HRV  metrics 
evaluated,  only  four  metrics  (DFA  Long,  FW,  SD1/SD2,  and  FD-L; 
three  of  which  were  calculated  by  non-linear  methods  to  assess 
RRI  irregularity  or  complexity)  supported  both  hypotheses. 


We  propose  that  a  metric  which  can  predict  low  tolerance 
in  an  individual  subject  prior  to  cardiovascular  collapse  would 
be  optimal  for  use  as  a  triage  assist  tool.  Identification  of  LT 
patients  is  particularly  important  as  the  first  responder  will  have 
less  time  to  initiate  effective  treatment  in  this  patient  population 
compared  with  HT  patients.  Upon  identifying  the  metrics  which 
supported  our  initial  hypotheses,  we  further  evaluated  DFA  Long, 
FW,  SD1/SD2,  and  FD-L  on  their  utility  as  triage  assist  tools  in 
individual  subjects.  First,  we  compared  the  individual  responses 
from  LT  subjects  at  presyncope  and  HT  subjects  at  submax  to 
determine  the  overlap  in  responses  between  the  two  groups.  Sec¬ 
ond,  we  assessed  the  accuracy  of  these  HRV  metrics  in  detecting 
LT  subjects  at  presyncope  from  the  total  subject  pool  (sensitivity 
and  positive  predictive  value). 

While  the  group  mean  values  of  DFA  Long,  FW,  SD1/SD2,  and 
FD-L  (Table  3)  distinguished  between  LT  at  presyncope  (unstable) 
and  HT  at  submax  (stable),  the  responses  in  individual  LT  and  HT 
subjects  varied  extensively  such  that  85-97%  of  the  LT  responses 
overlapped  with  those  of  HT  subjects  (Figures  3A-D).  In  addition, 
the  sensitivities  of  DFA  Long,  FW,  SD1/SD2,  and  FD-L  were  12, 6, 
24,  and  33%;  while  positive  predictive  values  were  50,  40,  67,  and 
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73%.  Considering  that  identifying  LT  subjects  by  chance  alone  has 
50%  sensitivity  and  50%  positive  predictive  value,  the  accuracy  of 
DFA  Long  and  FW  to  identify  LT  subjects  at  presyncope  was  not 
much  better  than  flipping  a  coin,  or  worse.  Essentially,  DFA  Long 
and  FW  did  not  adequately  identify  a  LT  individual  even  when  the 
subject  was  on  the  verge  of  cardiovascular  collapse. 

According  to  our  results,  FD-L  had  the  highest  positive  predic¬ 
tive  value  (73%),  and  thus  had  the  highest  potential  for  accurately 
assessing  hypovolemia  of  all  20  HRV  metrics  evaluated.  How¬ 
ever,  the  utility  of  FD-L  to  accurately  monitor  individual  trauma 
patients  may  be  limited  by  potentially  high  variability  in  individual 
responses.  In  the  present  study,  under  very  controlled  laboratory 
conditions,  the  sensitivity  of  FD-L  to  identify  a  LT  subject  was  only 
33%,  and  the  overlap  of  FD-L  measurements  from  LT  subjects  and 
HT  subjects  was  85%  (Figure  3D).  Overall,  these  results  indicate 
that  the  utility  of  FD-L  as  an  accurate  triage  assisting  tool  for  first 
responders  is  also  limited. 

Heart  rate  variability  has  been  studied  extensively  during  hem¬ 
orrhage  and  trauma.  HRV  decreases  with  hypovolemia  in  human 
LBNP  studies  (Cooke  and  Convertino,  2005;  Cooke  et  al.,  2008) 
and  in  animal  hemorrhage  experiments  (Batchinsky  et  al.,  2007b, 
2010).  HRV  has  also  been  studied  in  pre-hospital  environments 
during  transport  of  trauma  patients.  By  retrospective  analysis  of 
ECG  data  from  actual  trauma  patients  in  transport  to  hospital  care, 
a  number  of  investigators  (Cooke  et  al.,  2006a, b;  Batchinsky  et  al., 
2007a;  Cancio  et  al.,  2008;  Ong  et  al.,  2008;  King  et  al.,  2009)  have 
identified  several  HRV  metrics  that  are  associated  with  mortality 
or  the  need  for  a  LSI.  Importantly,  however,  all  of  these  analy¬ 
ses  were  based  on  evaluation  of  group  mean  data  only,  with  no 
consideration  of  the  appropriateness  of  application  to  individual 
patients. 

Despite  the  observation  of  depressed  HRV  in  subsets  of  trauma 
patients,  the  use  of  HRV  as  a  prognostic  tool  to  assess  the  severity 
of  hemorrhage  or  trauma  is  controversial  for  a  number  of  impor¬ 
tant  reasons.  Ryan  et  al.  (2010)  evaluated  the  ability  of  numerous 
HRV  metrics  to  track  hypovolemia  in  individual  subjects  under¬ 
going  simulated  hemorrhage  with  LBNR  They  observed  that  when 
group  means  were  evaluated  several  metrics  correlated  very  well 
with  reductions  in  stroke  volume  (r  >  0.87),  but  none  of  the  HRV 
metrics  consistently  correlated  with  changes  in  stroke  volume  in 
individual  subjects  (r  <  0.49). 

To  further  investigate  the  prognostic  relevance  of  HRV, 
Rickards  et  al.  (2010b)  evaluated  a  wide  range  of  HRV  metrics 
and  their  association  with  the  administration  of  LSIs  in  actual 
trauma  patients  who  all  had  normal  vital  signs  during  transport. 
These  patients  would  benefit  most  from  identification  of  an  early 
predictor  of  cardiovascular  collapse,  as  their  physiological  status 
could  not  be  accurately  determined  from  currently  available  stan¬ 
dard  vital  signs.  Of  the  HRV  metrics  studied,  only  one  metric, 
FD-L,  was  uniquely  associated  with  the  administration  of  a  LSI. 
However,  FD-L  variance  in  both  groups  of  patients  was  too  high 
to  accurately  determine  group  membership  (LSI  vs.  No-LSI)  on 
an  individual  basis,  and  the  number  of  false  negatives  identified 
with  FD-L  further  limited  the  power  of  this  metric  as  an  accurate 
indicator  of  LSIs  in  individual  trauma  patients.  Interestingly,  FD-L 
also  showed  the  highest  positive  predictive  value  for  identifying  LT 
subjects  in  the  current  study  but,  just  as  in  trauma  patients,  the 


overlap  between  individual  responses  from  LT  and  HT  subjects 
was  large  (85%),  and  we  found  the  sensitivity  (33%)  of  this  metric 
to  be  limited. 

There  are  several  technical  factors  associated  with  ECG  moni¬ 
toring  which  can  limit  the  usefulness  of  HRV  metrics  in  a  pre¬ 
hospital  setting.  First,  ECG  signals  must  display  normal  sinus 
rhythm  for  accurate  calculation  of  HRV,  but  trauma  patients  often 
develop  arrhythmias  such  as  premature  atrial  and  ventricular  con¬ 
tractions  (Sethuraman  et  al.,  2010).  Second,  ECG  signals  must 
have  a  low  level  of  random  noise,  which  usually  results  from 
motion  artifact,  an  unavoidable  consequence  of  patient  trans¬ 
port  in  the  pre-hospital  setting.  Third,  the  data  from  the  ECG 
signal  should  be  stationary,  i.e.,  a  relatively  stable  RR1  signal 
without  wide  fluctuations  from  interventions  or  patient  manipula¬ 
tion.  Unfortunately,  standard  pre-hospital  care  typically  requires 
extensive  patient  interventions  and  manipulation.  Other  factors 
which  can  also  contribute  to  ECG  non-stationarity  such  as  dis¬ 
ease,  age,  recreational  drugs,  medications,  alcohol,  smoking,  and 
postural  changes  are  widely  present  in  the  trauma  patient  popu¬ 
lation  (Bilchick  and  Berger,  2006).  Finally,  the  length  of  the  ECG 
data  set  required  for  accurate  measurements  varies  by  HRV  metric, 
and  can  range  from  100  to  800  RRIs;  which  can  be  an  inordinately 
long  time  in  a  trauma  patient  with  rapidly  changing  physiolog¬ 
ical  status  (Acharya  et  al.,  2006;  Rickards  et  al.,  2010a).  Thus, 
the  quality  of  HRV  measurements  are  optimal  when  they  can  be 
calculated  from  extended,  stable  ECG  recordings  under  standard¬ 
ized  and  very  stable  conditions;  these  conditions,  however,  are  not 
typically  encountered  when  first  responders  are  treating  trauma 
patients. 

Historically,  heart  rate  has  been  closely  monitored  in  trauma 
patients  because  of  the  widely  accepted  notion  that  the  magni¬ 
tude  of  tachycardia  reflects  the  degree  of  hypovolemia.  However, 
the  reliability  of  tachycardia  in  response  to  hypovolemia  has  been 
questioned  as  tachycardia  can  be  absent  in  many  trauma  patients 
despite  the  development  of  hypotension  associated  with  bleeding 
(Victorino  et  al.,  2003;  Brasel  et  al.,  2007).  In  addition,  the  results 
of  the  current  study  indicate  that  tachycardia  was  a  poor  indicator 
of  tolerance  to  hypovolemia  because  LT  subjects  at  presyncope  had 
lower  heart  rates  than  HT  subjects  at  greater  levels  of  central  hypo¬ 
volemia  (Table  3).  Based  on  the  clinical  doctrine  that  tachycardia 
signals  a  more  severe  state  of  hypovolemia  and  the  approach  of  car¬ 
diovascular  collapse,  the  HT  subjects  would  have  been  erroneously 
identified  as  more  hemodynamically  unstable  and  at  greater  risk  of 
developing  circulatory  shock  than  LT  subjects.  Although  a  change 
in  cardiac  rhythm  may  represent  an  adverse  clinical  status  during 
hemorrhage,  the  results  of  the  present  study  reinforce  the  unreli¬ 
ability  of  heart  rate  (and  subsequently  calculated  metrics  derived 
from  heart  rate  such  as  HRV)  for  patient  triage. 

In  summary,  the  results  of  the  current  study  are  consistent  with 
the  findings  that  ECG-derived  metrics  of  HRV  failed  to  provide 
reliable  information  about  clinical  status  in  individual  subjects 
during  progressive  reductions  in  central  blood  volume  similar  to 
those  experienced  during  hemorrhage  (Ryan  et  al.,  2010),  or  the 
need  for  LSIs  in  individual  trauma  patients  (Rickards  et  al.,  2010b). 
In  the  present  study,  heart  rate  and  HRV  metrics  derived  from  ECG 
signals  were  found  to  be  poor  indicators  of  LT  to  hypovolemia 
even  in  the  controlled  experimental  environment  of  simulated 
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hemorrhage.  This  study  raises  further  concern  that  monitoring 
heart  rate  or  calculated  derivatives  of  heart  rate  (e.g.,  HRV)  will 
not  reliably  identify  those  patients  who  are  least  tolerant  to  hypo¬ 
volemia  and  therefore  at  highest  risk  for  early  hemodynamic 
collapse  during  hemorrhage. 
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